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SUMMARY 

An equation was derived for the strain energy release rate, G, associated 
with local delamination growth from a matrix ply crack. The critical G c for 
edge delamination onset in [±25/90 2 ] s graphite epoxy laminates was measured and 
used in this equation to predict local delamination onset strains In 
[±25/90„] s , n = 4, 6, 8 laminates. A simple technique for predicting strain 
concentrations in the primary load bearing plies near local delaminations was 
developed. These strain concentrations were responsible for reduced laminate 
nominal failure strains in laminates containing local delaminations. The 
influence of edge delamination and matrix crack tip delamination on laminate 
stiffness and strength was compared. 
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area 

delamination size 

axial modulus of a laminate 

axial laminate modulus before delamination 

modulus of a locally delaminated cross section 

modulus of an edge delaminated laminate 

modulus of a matrix-ply-crack induced delaminated laminate 
strain energy release rate associated with delamination growth 
critical strain energy release rate for delamination onset 
local strain concentration 
gage length 

number of delaminations growing from a matrix crack 
number of plies 
applied axial load 

reciprocol of axial laminate modulus 
laminate thickness 

thickness of a locally delaminated cross section 
material volumne 
laminate width 
nominal axial strain 

nominal axial strain at onset of delamination 
in-si tu failure strain of primary load bearing plies 
nominal ultimate laminate strain at failure 
nominal axial stress 
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INTRODUCTION 


One commonly observed failure mode in laminated composite materials is 
delamination between the composite layers. Delaminations may form and grow 
under both static and cyclic tensile loading. The most common source of 
delamination is the laminate edge, where high interlaminar stresses are 
developed due to the mismatch in Poisson contraction of the individual plies. 
These edge delaminations typically occur between 90° plies and adjacent angle 
plies, with delaminations forming initially in a thumbnail shape, and rapidly 
becoming a delaminated strip that grows across the specimen width [1, 2], The 
strain energy release rate, G, associated with this edge delamination growth, 
and the contributions of interlaminar tension and shear to G, have been 
determined previously for a variety of layups [1-4]. The influence of edge 
delamination on tensile stiffness and strength has been determined for a 
variety of layups as well [1, 5], 

Another source of delamination is matrix ply cracks running parallel to 
the fibers in a ply. The interlaminar stresses that develop in the ply inter- 
face at matrix crack tips may cause local delaminations to form and grow [3, 

6, 7]. The purpose of this investigation was to determine delamination onset 
strains for these local de laminations, and to determine their influence on 
tensile stiffness and strength. To this end, a strain energy release rate 
analysis of local delaminations was developed. Then, a simple technique for 
calculating strain concentration factors in the primary load bearing plies 
near localized delamination was developed. Finally, the influence of edge 
delamination and matrix crack tip delamination on laminate stiffness and 
strength were compared. 
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STRAIN ENERGY RELEASE RATE 


Analysis 

For an elastic body containing a crack that grows under a constant 
applied load, P, the strain energy release rate, G, is given by [8] 

P 2 dC 

g - — *k ( i ) 


where C is the compliance, i.e., C — 6/P, and A is the crack surface area 
created. A similar expression may be written for G in terms of the remote 
stress, a, and the reciprocal of modulus, s, where S = e/a, by substituting 



into equation (1). This yields 


G 



2 dS 
dA 


(3) 


where dS/dA is the rate of change in S as the flaw extends, and V is the 
volume of the body, i.e., V = wt£ . 

In this study, the elastic body is a composite laminate containing a 
matrix ply crack through the thickness of n off-axis plies, with delami- 
nations forming at the matrix crack tip and growing in the ply interfaces 
(Fig. 1). For simplicity, the strain energy release rate associated with the 
growth of delaminations from a single matrix crack will be considered. 

In order to evaluate dS/dA in equation (3), and hence G, an equation for 
laminate compliance as a function of delamination size was developed. 

Figure 2 illustrates a composite laminate containing delaminations 
growing from a matrix ply crack. The composite gage length, £, is divided 
into a locally delaminated region, a, and a laminated region, £ - a. 
Assuming the composite displacements are the sum of the displacements in these 
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two regions, and the total load, P, is equal to the loads carried by the two 
regions individually, then using Hooke's law, AJ l = Pi/AE, yields 

S = 1 = Alam U ~ a) a 

E * . ^LAM^LAM a ld e ld 

where A LAM and Ej^ are the cross sectional area and modulus of the 
laminated region, and A^ and E^ are the cross sectional area and modulus 
of the locally delaminated region. Each of the areas in equation (4) repre- 
sent only the cross sectional area that carries the applied load, hence. 



A 


LAM 


wt 


(5) 


A 


LD 


wt 


LD 


( 6 ) 


where w and t are the laminate width and thickness, respectively , and t ■ 
is the thickness of the locally delaminated region that carries load, i.e # , 
the thickness of the uncracked plies. Substituting equations (5) and (6) into 
equation (4) yields 



Hence, the laminate compliance is a linear function of delamination size, a. 

Returning to Figure 2, the strain energy release rate associated with the 
growth of delaminations from a matrix ply crack can be calculated by assuming 


V = tw i 
A = mwa 
dA = mw da 


(8) 


where m is the number of delaminations growing from the matrix ply crack. 
For the case illustrated in Figure 2, m = 2, but for a delamination growing 
from a cracked surface ply, m = 1. Substituting equations (8) into equa- 
tion (3) and differentiating equation (7) yields 
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G = 


-_\ 

2m \ t E tE / 
\ LD LD ^LAM / 


( 9 ) 


or similarly, in terms of the applied load p, 
2 


G = 


( t LD E LD ' “LAM ) 


( 10 ) 


2mw \ LD LD “"LAM 
Hence, as indicated in equations (9) and (10), the strain energy release rate 
is independent of delamination size. The magnitude of G depends only on the 
laminate layup and thickness, the location of the cracked ply and subsequent 
delaminations (which determines E^, t^, and m), the applied load, P, and 

the laminate width, w. 


Experiments 

Previously, Crossman, Wang, and Law [2] performed tests on T300/934 
[±25/90 n ] g laminates, where n varied from 1/2 to 8. They recorded the load 
at the onset of transverse cracking in 90° plies, at the onset of delamina- 
tion, and at failure. They used dye penetrant enhanced X-radiography to 
confirm the onset and extent of damage. Their results are illustrated in 
Figure 3, which shows the tensile strain at the onset of transverse cracking 
in the 90° plies, at the onset of delamination, and at final failure, as a 
function of the number of 90° plies, n, in the [±25/90 n ) g laminates. 

A significant reduction occurs in both the delamination onset strain and 
failure strains for the n > 4 laminates. Figure 4 shows a schematic of 
fracture sequences across the laminate width and at the edge (a) just before 
the onset of delamination, (b) just after the onset of delamination, and 
(c) just before final failure. The sequence illustrated in column (b), as 
n increases from 1/2 to 8, indicates a transition between n = 3 and 4 from 
thumbnail-shaped edge de laminations to local de laminations growing from a 
90 matrix ply crack. Hence, these data are useful for checking the strain 
energy release rate analysis for local de laminations. 
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Delamination Onset Prediction 


First, a critical value of G c was calculated from a previously derived 
equation for edge delamination [1] f 

0 E *) do 

using the measured delamination onset strain for the [±25/90 2 ] g laminate. The 
totally delaminated modulus, E*, was calculated using laminated plate theory 
and the rule of mixtures [1] 

q 

e. t. 




where E* = modulus of the delaminated laminate 

q = number of sublaminates formed by the delamination 
E i = modulus of the ith sublaminate formed by delamination, and 
= thickness of the ith sublaminate 

As shown in the edge views of Figure 5(a), edge delaminations typically form 
at interfaces between 90° plies and adjacent angle plies, shifting from one 
interface to its symmetric counterpart along the gage length. Previous work 
[1,5] has indicated that measured stiffness loss agreed well with E* calcu- 
lations from equation (12) when delaminations were modeled at both interfaces 
between angle plies and 90° plies and the sublaminated stiffnesses, E^ were 
calculated from laminated plate theory. For the t±25/90 n ) s family of layups. 


with delaminations modeled in the -25/90 interfaces, 


E* = 


4 E (±25) + 2n E 90 

s 

(4 + 2n) 


( 13 ) 


which, for the n = 2 laminate used to measure G , becomes 

c 
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E* = 


E (±25) + E 90 

s 


(14) 


Using lamina properties from [9] + 


E 11 = 

134 GPa 

(19.5 

Msi ) 

E 22 = 

10.2 

GPa 

(1 .48 

Msi ) 

G 1 2 = 

5.52 

GPa 

(0.8 1 

Msi) 
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in equations 

(14) and (11) yielded a 

Then, equation (11) 

i ■■ 

was 

solved 

for e 


e = 
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\J 


2 G 


t(E - E*) 
LAM ' 


(15) 


and this G c was used to predict £ c for edge delamination onset in the 

other [±25/90 n ) g laminates. As shown in Figure 6, good agreement was observed 

for n = 1 and 3 where edge delamination was observed (Fig. 4), but the 

predicted e was much higher than measured onset strains for the n > 4 
c 

laminates where delaminations formed at 90° ply cracks. 

A similar procedure was used to predict delamination onset from 90° 

matrix ply cracks. Because no significant stiffness loss was observed before 

local delamination onset [2], Hooke's law a = E e was substituted into 

LAM 

equation (9), and equation (9) was solved for e yielding 


e = 
c 


2mG 


E t 
LAM 


i 


l 


E LD fc LD E LAM t , 


( 16 ) 


tThese lamina properties were chosen because they resulted in the best 
comparison of values calculated from laminated plate theory and 

measured modulus reported in reference 2. 
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Then, the G Q from the n = 2 laminate was used in equation (16) to predict 
e c for local delamination onset in the other [±25/90 n l g laminates. For these 
local delaminations growing from a single matrix crack through the thickness 
of the entire set of n 90° plies, t^ D and E^p were simply the thickness 
and modulus, respectively, of the four remaining ±25° plies, and m was equal 
to 2 (Fig. 2). As shown in Figure 6, the delamination onset strains agreed 
fairly well for the n > 4 laminates, and captured the trend of decreasing 
with increasing thickness of the cracked 90° plies. 

The results in Figure 6 indicate that a strain energy release rate 
analysis will account for delamination onset for both sources of delamination, 
as illustrated in Figure 6, and account for the observed thickness dependence. 
However, the simple equations for total energy release rate shown in Figure 6 
may not be sufficient for correlating all possible combinations of layups and 
delamination locations. Previous work [4, 10, 11] has indicated that for 
brittle resin composites with relatively low values of G c , similar to values 
measured in this study for T300/934, the mode X component of G, due to inter- 
laminar tension only, controls the onset of delamination. This mode I depen- 
dence may explain why the transition from edge delamination to matrix ply 
crack induced delamination, that could be inferred from the predictions in 
Figure 6 to occur at n = 2, does not occur until n = 4. To verify the 
fracture mode dependence, the boundary value problem for local delamination 
must be formulated and solved to determine the contributions of interlaminar 
tension and shear to the total G calculated in equations (10) and (11). 

TENSILE STIFFNESS AND STRENGTH 

As shown previously (Fig. 3), local delaminations that grew from matrix 
ply cracks in [±25/90 n l s laminates (n > 4) formed at lower strains than edge 
delaminations (n < 3). Figure 3 also shows that the laminate failure strains 
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for [±25/90 n ] s laminates that experienced edge delamination (n < 3) are nearly 
identical to the failure strains of [±25 ] 2s laminates (n = O'). But the 
laminates that experience local, matrix-ply crack induced delamination 
(n > 4) failed at lower strain levels. The nominal ultimate failure strain, 
e nu * decreased as n increased. Furthermore, these n > 4 laminates failed 
soon after the onset of local delamination. This data indicates that local 
delaminations from 90° ply cracks create local strain concentrations in the 
primary load bearing ±25° plies that precipitates early laminate failures. 
Therefore, an analysis was developed to determine local strain concentrations 
in the load bearing plies in the vicinity of local delaminations. 


Local Strain Concentration Analysis 

Figure 7 shows through- thickness free-body diagrams for [±25/90 n J g 
laminates where (a) edge delamination has been modeled in the -25/90 inter- 
faces, (b) where local delamination has been modeled in the -25/90 interfaces 
growing from a 90^ matrix ply crack, and (c) where both edge and local 
delamination are modeled. Assuming the load at the top (T) and bottom (B) 
cross section (as shown in Fig. 7) must be equal, then substituting Hooke's 
law, yields an expression for the local strain concentration, K, from the 
top to the bottom through- thickness cross section. 


K = 


e e t 

B _ T T 

e e t 

T B B 


(17) 


For case (a) where only edge delamination occurs, the modulus of the top and 
bottom sections are identical, and are equal to E* as calculated by 
equation (13). Similarly, because all the plies continue to carry load, the 
thickness of the top and bottom cross sections are identical. Hence, edge 
delamination does not result in local strain concentrations in the individual 
plies. This conclusion is consistant with the observation that the n < 3 
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laminates that edge-deldminate had the same failure strains as t±25]2 s 
laminates. However, for case (b) where delaminations extend from a 90° ply 
crack, local strain concentrations do occur. The modulus and thickness of the 
top cross section are those of the original laminate. The modulus of the 
bottom cross section is simply the modulus of a [±25] s laminate. The thick- 
ness of the bottom cross section is simply the thickness of the four ±25° 
plies, because the isolated (cracked and delaminated) 90° plies no longer 
carry load. Hence, the isolated 90° plies are modeled as if they have been 
removed from the laminate (Fig. 5(b)). Table 1 shows the resulting local 
strain concentrations, K, in the load bearing ±25 plies calculated for the 
n = 4, 6, 8 laminates. These local strain concentrations should account for 
the reduction in [±25/90 n ] s laminate nominal failure strains observed as n 
increases from n < 3, where only edge delamination is observed, to n > 4, 
where local delamination is observed. However, as shown in Figure 4, the 
n > 4 laminates exhibit a combination of edge and local delamination before 
final fracture, even though the local delamination occurs first. Hence, the 
free body in Figure 7(c) may be more appropriate for evaluating local strain 
concentrations. For this case (c), the modulus of the top cross section is 
equal to E* from equation (13) and the thickness is the original laminate 
thickness. The modulus and thickness of the bottom cross section is simply 
the modulus and thickness of a [±25] g laminate. The resulting local strain 
concentrations for this case (c) are also shown in table 1 . 

These local strain concentrations were used to predict the failure strain 
of the n = 4, 6, 8 laminates. Because the n < 3 laminates contained only 
edge delamination at failure, which does not result in local strain concentra- 
tions, their laminate nominal failure strains were also the in-situ failure 

strain , e , of the ±25° plies. Hence, a mean laminate failure strain of 
F 
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0.0065 was assumed to be the failure strain of the primary load bearing ±25° 
plies (Fig. 3). Then, the n = 4, 6, 8 laminate failure strains were pre- 
dicted as e = e /K, where K was the local strain concentrations listed in 
nu r 

table 1(c). The measured and predicted failure strains agreed well (Fig. 8). 
The trend of decreasing nominal failure strain with increasing thickness of 
the cracked ply is evident. 

Stiffness Loss 

l 

Because local delaminations growing from matrix ply cracks immediately 
result in local strain concentrations when they form, laminate failure may 
occur soon after the formation of local delaminations. This was indeed the 
case for n > 4 laminates, as shown in Figure 3. Laminate failure occurred 
at, or just above, the onset strains for local delamination. If, however, 
these local delaminations were able to grow over a significant portion of the 
gage length, an associated stiffness loss would be observed. Substituting 
a = l into equation (7) yields 



which represents the laminate modulus if the delaminations that formed at a 
matrix ply crack grew over the entire gage length. This totally delaminated 
modulus, E* , is not the same as E* calculated for edge delamination from 
equation (13), because the 90° plies become isolated (cracked and delaminated 
on both sides as shown in Fig. 5(b)) for delaminations growing from matrix ply 
cracks. For example. Figure 9 shows the modulus (normalized by Ej^^) of 

f±25/90 n ) g laminates (n = 1/2 to 8) completed delaminated at the -25/90 * 

interfaces. For edge delamination, the largest modulus decrease is 23 percent, 
and occurs at n = 2 . However, if the delaminations grew from a 90® matrix 
ply crack, resulting in isolation of the 90® plies, then the laminate modulus 
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decreases monotonically with increasing n. For n = 8, a 49 percent decrease 
in laminate modulus would occur. 

Comparison of Edge and Local Delamination 

Figures 10 and 11 summarize the differences between edge delamination and 
matrix crack tip delamination. As shown in Figure 10, edge delaminations form 
at a constant strain £, governed by equation (11). As edge delaminations 
grow across the width, the load carried by the delaminated strips near the 
edges is reduced. Eventually, when delamination has grown through the width, 
the entire laminate carries a reduced load at the same strain, and has a 
modulus E* determined by equation (12). Hence, edge delamination is 
inherently a stable fracture process. In contrast. Figure 11 shows that local 
delaminations form at matrix crack tips at a constant load P, governed by 
equation (10). As local delaminations begin to grow from the matrix crack tip 
and isolate the cracked plies, the strain in the uncracked plies near the 
locally delaminated region increases, with local strain concentrations estimated 
by equation (17). These local strain concentrations may result in failure of 
the uncracked plies that carry the load before extensive delamination growth 
occurs. If these delaminations were able to grow over the entire gage length, 
the uncracked plies would have an increased uniform nominal strain, and the 
laminate would have a reduced modulus determined by equation (18). Hence, 
delamination growth from matrix ply cracking may be unstable, leading to 
immediate laminate failure as was observed for t±25/90 n ) s n = 4, 6, 8 
laminates; but, in other laminates, local delamination may result in some 
stable growth before failure, or before other interfaces delaminate, which may 
further increase local strains and cause failure. This latter behavior has 
been observed under cyclic tensile loading of quasi-isotropic laminates [12]. 
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SUMMARY AND CONCLUSIONS 


A simple equation for the strain energy release rate, G, associated with 
local delaminations growing from matrix ply cracks was developed. This 
equation was then used to predict the local delamination onset strains in 
[±25/90 n ] g laminates where delaminations grew from matrix cracks in the 
90° plies. This analysis correctly predicted the trend of lower delamination 
onset strains with increasing thickness of the cracked ply. 

A simple technique was developed for calculating local strain concentra- 
tions in the primary load bearing plies near localized de laminations. These 

strain concentrations were used to successfully predict [±25/90 ] laminate 

n s 

failure strains. 

Based on these observations, and previous work, the following conclusions 
have been reached: 

(1) Strain energy release rate, G, is a reasonable generic parameter for 
characterizing the interlaminar fracture behavior of composites. This 

G parameter accounts for the observed thickness dependence for both edge 
delamination and local delamination onset, and allows correlation between 
these two types of delaminations. 

(2) Edge delamination is a stable fracture process that may reduce laminate 
modulus, and hence influence tensile strength, but will not cause 
premature laminate failures. 

(3) Local delaminations growing from matrix cracks, however, create local 
strain concentrations that may lead to nominal laminate failure strains 
below the in-situ failure strain of the primary load bearing plies. 
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TABLE 1 : LOCAL STRAIN CONCENTRATIONS IN ±25° PLIES OF 

t±25/90 n ] g LAMINATES 

K - ^ptip/Egtg 


n 

(a) Edge Delam Only 

(b) Local Delam Only 

(c) Combined 

4 

1 .00 

1 .66 

1 .29 

6 

1 .00 

1 .81 

1 .43 

8 

1 .00 

1 .95 

1 .57 
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EDGE VIEW FRONT VIEW 


Fig. 1 — Schematic of local delaminations growing from matrix cracks. 






[±25/90 n ] T300/934 



Fig. 4 — Schematic of fracture sequences in [ ±25 /90 n ] ^ laminates. 
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Fig. 5 — Sources of delamination in unnotched laminates. 
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Fig. 7— Free body diagrams of local through thickness sections of [±25/90 n ] laminates 
containing delamination. s 
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Fig. 8— Prediction of [±25/90 n ] g laminate failure strains. 







1. Report No. USAAVSCOM* !*• Government Accession No. 

NASA TM-85728, TR-83-B-6 


4. Title and Subtitle 

ANALYSIS OF LOCAL DELAMINATIONS AND THEIR INFLUENCE 
ON COMPOSITE LAMINATE BEHAVIOR 


7. Author(s) 

T. Kevin O'Brien 


9. Performing Organization Name and Address 

Structures Laboratory 

USAAVSCOM Research and Technology Laboratories 
NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

and 

U.S. Army Aviation and Systems Command 
St. Louis, MO 63166 

15. Supplementary Notes ™ " 

Presented at the ASTM Symposium on Delamination and 
Pittsburgh, Pennsylvania, November 8-10, 1983. 

16. Abstract 


3. Recipient's Catalog No. 


5. Report Date 

January 1984 


6. Performing Organization Code 

506-53-23-05 


8. Performing Organization Report No. 


10. Work Unit No. 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 

14. Army Project No. 

1L161102AH45 


Debonding of Materials, 


An equation was derived for the strain energy release rate, G, associated 
with local delamination growth from a matrix ply crack. The critical Gj. for 
edge delamination onset in [±25/902 ] s graphite epoxy laminates was measured and 
used in this equation to predict local delamination onset strains in [±25/90n]s> 
n = 4, 6, 8 laminates. A simple technique for predicting strain concentrations 
in the primary load bearing plies near local delaminations was developed. These 
strain concentrations were responsible for reduced laminate nominal failure 
strains in laminates containing local delaminations. The influence of edge 
delamination and matrix crack tip delamination on laminate stiffness and strength 
was compared. 


17. Key Words (Suggested by Author(s) ) 

Graphite epoxy materials Fatigue 

Delamination Composite 

Matrix ply cracks materials 

Strain concentration Interlaminai 

Strain energy release rate fracture 


19. Security Classif. (of this report} 

Unclassified 


18. Distribution Statement 


Unclassified - Unlimited 


Subject Category 24 


20. Security Classif. (of this page) 

Unclassified 


21. No. of Pages 

29 


22. Price* 

A03 


For sale by the National Technical Information Service, Springfield, Virginia 22161 









f 



